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Abstract: Unified Dark Matter models describe Dark Matter and Dark Energy as a 
single entity which is, in the simplest case, embodied in a perfect barotropic fluid. It is 
a well-established fact that small adiabatic perturbations of Unified Dark Matter have an 
evolution characterised by oscillations and decay which provide predictions on the Cosmic 
Background Radiation anisotropies which are in poor agreement with observation. In this 
paper we investigate the generalised Chaplygin gas and we find that the Integrated Sachs- 
Wolfe effect excludes the model for 10~^ < a < 350. We discuss the implications on the 
background evolution of the Universe if large values of a are considered. In this case, the 
Universe expansion mimics a matter-dominated phase abruptly followed by a de Sitter one 
at the transition redshift ztr- Thanks to an analysis of the type la supernovae Constitution 
set we are able to place ztr = 0.22. 

Keywords: generalised Chaplygin gas. Integrated Sachs- Wolfe effect, matter-radiation 
equivalence scale, Meszaros effect, SNIa Constitution set . 
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1. Introduction 

The generalised Chaplygin gas (hereafter gCg) is a cosmological model which has raised 
some interest thanks to its ability to describe alone the dynamics of the entire dark sector 
of the Universe, i.e. Dark Matter and Dark Energy. The gCg is a perfect fluid characterised 
by the following equation of state: 

P = -— , 1-1 
P 

where p is the pressure, p is the density and A and a are positive parameters. 

In the setting of the Friedmann-Lemaitre-Robertson- Walker cosmological theory, one 
easily integrates the energy conservation equation and finds that the gCg density evolves 
as a function of the redshift z as follows: 

1 



P 



(1.2) 



where i? is a positive integration constant. Equation ([T]^) interpolates between a past 

[z ^ 1) dust-like phase of evolution, i.e. p ~ (1 + z)^, and a recent {z < 1) one in which p 

1 

tends asymptotically to a cosmological constant, i.e. p ~ A'^+°' . 

As cosmological model, the gCg was first investigated in because it had already 
raised some attention, e.g. in connection with string theory. In particular, the equation of 
state (|ri|) can be extracted from the Nambu-Goto action for d-branes moving in a (d + 2)- 
dimensional spacetime in the lightcone parametrisation, see Q. Also, the gCg is the only 
fluid which, up to now, admits a super symmetric generalisation, see M, H. 
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From the point of view of cosmology, a possible unification picture between Dark Mat- 
ter and Dark Energy is particularly appealing, especially in connection with the so-called 
cosmic coincidence problem (see H] for much detail about the latter). This motivation 
prompted an intensive study of the gCg and of those models that have the property of 
unifying (dynamically) Dark Matter and Dark Energy. Such models are called Unified 
Dark Matter (UDM). 

It must be pointed out that in a conventional cosmological model (i.e. a model in which 
Dark Matter and Dark Energy are separated entities) Dark Matter has the primary role 
of providing for structure formation whereas Dark Energy has to account for the recently 
observed accelerated expansion of the Universe [P, |7|. The fundamental task for the gCg 
is then to play both these roles. 

To understand if this is the case, the gCg parameters space {A, a) has been intensively 
analysed in relation with observation of Large Scale Structures (LSS), Cosmic Microwave 
Background (CMB), type la Supernovae (SNIa), X-ray cluster gas fraction and Baryon 
Acoustic Oscillations (BAO), see for example §, |, |ll], ||, ||, 0, |l8|. 



The most instructive constraints are about the parameter a and come from the LSS 
matter power spectrum and the CMB angular power spectrum analysis. In the former, the 
resulting best fit value is a < 10^^, |13|. One can think that this narrow constraint is a 
flaw of the gCg model because if we take the limit a — )• in the equation of state Eq. (p. . Ij) , 
then we reproduce exactly the ACDM model (if A assumes the value of the cosmological 
constant energy density). 



The authors of [|19[ have found a loophole in this degeneracy problem. Indeed, they 
have shown that if we consider a baryon component added to the gCg and we compute 
the matter power spectrum for the baryonic part alone, it turns out that the latter is only 
poorly affected (at the linear perturbative regime) by the presence of the gCg. The model 
constituted by gCg plus baryons is therefore in good agreement with LSS observation for 
all a e (0,1). 

On the other hand, even if the "baryon loophole" allows to circumvent the very narrow 
constraint on a coming from LSS analysis, it cannot prevent the tight one coming from 



CMB analysis and due in particular to the Integrated Sachs- Wolfe (ISW) effect, [|14|, |15 



|16| , 17, 18, In more detail, if we take into account an ordinary Cold Dark Matter 
(CDM) component and a baryonic one together with the gCg, we flnd that a < 0.2, 
|T^. Withdrawing the CDM component, we lower the bound by an order of magnitude: 
a < 10~^, [|l^]. Finally, for the case of the pure gCg we flnd an even tighter constraint: 
a < 10-^, H. 

Taking into account all these results, we may conclude that the gCg model is viable 
only when it is very similar, if not degenerate, to the ACDM. However, it must be pointed 
out that in the major part of the literature on the subject, the preliminary constraint q < 1 
is assumed. The reason for this assumption resides in the form of the gCg adiabatic speed 
of sound, which is the following: 

- ? = ^^^^TTTT • (1-3) 

dp ^ + 5(l + ^)3(-+l) 
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If z — )• —1 then Cg — ?■ a. Therefore, causaUty would require a < 1. However, the causality 
problem for UDM models should rather be addressed in terms of a microscopic theory, see 
|21]. In the particular case of the gCg, for a > 1 the authors of [^] examine in some detail 
the causality issues and develop a suitable microscopic theory in which the signal velocity 
never exceeds the speed of light. 

Given this, it is natural to wonder how the gCg model behaves in the "forbidden" range 
a > 1. In the literature few papers carry out such analysis, for example 23, 24, 

An interesting result is for instance that, in the framework of the "baryon loophole" above 
mentioned, i.e. considering the model constituted by the gCg plus baryons, the agreement 
of the baryon power spectrum with observation increases for a > 3, see [^]. In [23| and [24| 
the authors have recently confirmed this behaviour. On the other hand, CMB constraints 
do not change significantly from the case a < 1, in the sense that the parameter is again 
constrained to be very small, see |17| and [25|. 

Even in those papers which analyse the gCg for a > 1, the parameter space is probed 
up to a maximum finite value. For example, q < 6 in [^] or a < 10 in [^]. Since the 
a > 1 case has proved to be promising, at least in relation with LSS analysis, in this 
paper we dedicate ourselves to investigate the extreme limit of the gCg, i.e. the behaviour 
of the model for very large values of a. Our analysis is principally based on the ISW 
effect because, as we have discussed in the above and as the authors of |2^ have shown, it 
provides the strongest constraints for a UDM model. 

To this purpose, in section II we briefly outline the basic equations describing the ISW 
effect and, inspired by |2C], we present a simple method based on the Meszaros effect which 
we employ in Section III to find a qualitative constraint for large values of a. Indeed, in 
section III we analyse the Jeans wavenumber of gCg perturbations and we find that if 
a > 350 the ISW effect can be potentially small. In Section IV we then confirm the results 
of Sec. Ill by directly calculating the ISW effect in the limit a — t- oo and showing that it is 
smaller than the one for a = 0, i.e. for the ACDM (the calculation of the ISW effect for the 
ACDM has been performed for the first time by Kofman and Starobinsky in 1985, see [^]). 
In section V we address the behaviour of the background expansion for large values of a 
and find that the evolution is characterised by an early dust-like phase abruptly interrupted 
by a de Sitter (dS) one. We then analyse the 157 nearby SNIa of the Constitution set and 
find that the transition between the two phases takes place about a redshift ztr = 0.22, 
which is much more recent than the transition at ztr = 0.79 to the accelerated phase of 
expansion in the ACDM model. The last section is devoted to discussion and conclusions. 



2. The ISW effect and the constraining method 

In this paper we discuss only adiabatic perturbations and we assume a flat spatial geometry. 
The ISW effect contribution to the CMB angular power spectrum is given by the following 



formula |27]: 



2^ + lwsw_ 1 /■°°dfc,,|GKr/o,fc)|' 



47r ^' - 2^ I 2Z + 1 ' ^^-^^ 
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where / is the multipole moment, k is the wavelength number, r/o is the present conformal 
time and 

^ = 2 y (r/, k) ii [k im - r,)] dv (2.2) 

is the fractional temperature perturbation generated by the ISW effect, where 77* is the last 
scattering conformal time, ^ {rj, k) is the Fourier transformed gravitational potential and 
j; is the spherical Bessel function of order I. The prime denotes derivation with respect to 
the conformal time rj. For a more detailed description of the perturbations equations and 
of the integrals (^^1-2. 2) we refer the reader to |27, 29, pO[ . 



Following [^9[, consider the Fourier transformed evolution equation for the gravita- 
tional potential: 

9" 

u" + k^clu - — u = , (2.3) 



where: 



and Cs, p and p are, respectively, the adiabatic speed of sound, the energy density and 
the pressure defined for a generic cosmological model. We have chosen here units such 
that SttG = c = 1. For simplicity, here and in the following we do not explicit the (r?,A;) 
dependence for u and <^ and the rj dependence for c^, 6*, p, p and z. They will therefore be 
always implied, unless otherwise stated. 
In Eq. (|2]|) define 

as the square Jeans wavenumber. In general, when wavelengths smaller than the Jeans scale 
enter the Hubble horizon, they start to oscillate, affecting the CMB power spectrum and 
the matter one in ways not compatible with observation. For UDM models in particular, 
the authors of ||2^ have found a contribution to the ISW effect proportional to the fourth 
power of the speed of sound which generates in the CMB angular power spectrum a growth 
proportional to until / ~ 25 (the value / ~ 25 is related to the equivalence wave number 
keq which we discuss later). See [^] and also for more detail. 
Consequently, if we take into account only those scales for which 

k^ <kj, (2.6) 

then the gravitational potential does not oscillate. 

Of course if Cg = 0, as for the ACDM model, the Jeans scale vanishes and condition 
( |2.6| ) is satisfied for every k at any time (remember that fcj is time dependent). On the 
other hand, this is not the only possible scenario. In fact, the cosmological scales which 
are important for the CMB and structure formation are those which entered the Hubble 
horizon after the matter-radiation equivalence epoch. Those which entered the horizon 
earlier had been damped by the dominating presence of radiation. This effect is known as 



Meszaros effect g |3|, |3| 
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If we require that the relevant scales, i.e. k < kcq, must satisfy condition ( ^.6| ) we 
obtain: 

k'e, < 1^ • (2.7) 



C; 



This relation can be used to infer qualitative constraints upon a generic cosmological model. 
In the next section we will make use of it to find constraints on a. As a remark, a scenario 



for which Eq. (2^) holds true without demanding a vanishing speed of sound is the so-called 
fast transition, introduced and investigated in great detail in [03]. 



3. Constraints on the generalised Chaplygin gas 

The equivalence wavenumber feeq has the following form: 

tt2 

,2 _ eg ,x 

(1 + Zeq) 

where Heq is the Hubble parameter evaluated at the equivalence redshift Zeq- From the 5- 
year WMAP observation^ , the best fit values are Zeq = 3176tJ|J and keq = 0.00968±0. 00046 
h Mpc~^, where h is the Hubble constant in 100 km s~^ Mpc~^ units. See also [35, 



The Hubble parameter is related to the Universe energy content by Friedmann equa- 
tion, which for the pure gCg model has the following form: 

U = [i + (1 - i) (1 + zf ^ , (3.2) 
Hq L J 

where A = A/ (A + B) and Hq is the Hubble constant. Notice that A = —wq, where wq is 
the present time equation of state parameter of the Universe. 

Let ztr be the redshift at which the accelerated phase of expansion begins. From 
Eq. ( |3.2| ), we calculate the following relation between A and ztr^ 

i = ^i±^^I^^. (3.3) 

2+(l + Ztr)'^"+l^ 

From now on we make use of {ztr, a) as independent parameters. 

Plugging Eqs. 1^), (jT^), and (jsj) into the definition (|2j), we plot in 

Fig. m the Jeans wavenumber as function of a for fixed z = and ztr = 0.79 and the 
equivalence wavenumber, computed from Eq. (|3.2| ), as function of a for a fixed Zeq = 3176. 



The value we have chosen for the transition redshift Ztr is the WMAP5 best fit, see [^, |36[. 

The most intriguing feature of Fig. || is that the Jeans wavenumber has a minimum 
value for a « 1. For sufficiently small or sufficiently large values of a it grows and equals 
fceq- From our numerical computation we have inferred the following constraints: a < 10^^ 
and a > 250. 



^http : //lambda. gsfc .nasa. gov/ 
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Figure 1: fcj (black curve) and fcoq (red "quasi-horizontal" line) as functions of a. kj is evaluated 
at z = and ztr = 0.79, while fcoq is evaluated at Zoq = 3176. The wavenumbers are in units h 
Mpc"^ 



It is also possible to obtain these results from approximations, but analytically. Indeed, 



when a <C 1 we expand in Taylor series /Cgq and akj and find: 



•-eq 



ak'j 



{l + z, 



4 + (1 + Ztr 



2 + (1 + Ztr) 
3 



+ 0{a) , 



0{a) 



(3.4) 



(3.5) 



4(1 + Ztr)' 

To leading order in a, equating the above expressions gives the upper bound for small 
values of a. For Zeq = 3176 and ztr = 0.79: a < 10~^. 



Now expand kg„ in Taylor series for a » 1: 



"'eq 



l + .eq ^^(1 



(3.6) 

(1 + ^tr)^ ■ > " ' 

The corresponding expansion for kj is less immediate. Making use of the asymptotic forms 
of the speed of sound and of the Hubble parameter, which we will give in Eq. ( [4.2D and in 
Eq. (^), we find: 

6 



ki 



where we have defined 



,3a 



9a 



« (1 + Ztr 
1 



+ 0^ 



for X > 1 



(3.7) 



4 (l + z) 



2 [l + O (x^")] for X < 1 



X 



l + z 

1 + Ztr 



(3.8) 
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From Eq. (|3.7| ), to leading order in 1/a the Jeans wavenumber is an exponential function 
of a for X > 1; for x < 1, the expansion can be performed only with respect to x'^" and, to 
leading order, k'j grows linearly with a. We take into account the latter instance, equate 
( |3.7| ) to ( |3.6| ) (each expression considered to its leading order) and find a > 250 (for z = 0, 
= 3176 and Ztr = 0.79). 

Note that we have neglected the CDM, the baryon and the radiation components. 
Neglecting the CDM component is reasonable, because the gCg model aims to an unified 
description of Dark Matter and Dark Energy. Adding a CDM component would then spoil 
its purpose. Moreover, neglecting radiation is also reasonable, because the minimum value 
of the Jeans wavenumber is at late times {z w 0), where radiation is subdominant. For 
what concerns baryons, their presence would have the effect of lessening the average speed 
of sound, increasing thus the Jeans wavenumber and smoothing the constraints we have 
found. Nonetheless, at late times the baryon component is also subdominant with respect 
to the gCg one, so we can reasonably neglect it. 

On the other hand, in the calculation of k'^^ we are not allowed to neglect both radi- 
ation and baryons, because at Zgq = 3176 they are important. Therefore, since at early 
times the gCg and the ACDM model are indistinguishable, we use the WMAP5 result 
^eq = 0.00968 h Mpc~^ and we find that a < 10~^ and a > 350. As we expected, taking 
into account the radiation contribution has the effect of increasing the lower bound for 
large values of a. 



4. Calculation of the ISW effect for the extreme limit of the generalised 
C haply gin gas 



In this section we calculate the integrals {2A_) and ( |2.2D in the limit a — t- oo. To this 



purpose, in place of Eq. (2^), we employ the evolution equation for the gravitational 
potential namely 

+ m{i + c^) + {in' + n'^ + m^cl + k^c^) $ = o , (4.1) 

where Cg is the gCg adiabatic speed of sound defined in (|1.3D and % = a' /a is the Hubble 
parameter written in the conformal time (see [^] for more detail). In the limit of very 
large values of a, from Eq. ( |1.3|) together with Eq. ( |3.3|) , the square speed of sound has the 
following asymptotic behaviour 



r a 



[x-3° + O (x-6")] for X > 1 



(4.2) 



^ a [1 - 2x3° + O (x6°)] for X < 1 
with Cg = q;/3 for x = 1. Friedmann equation (3.2) becomes: 

ln2 _ / 1 



El 



1 + 



+ 



a 



a" 



. 2x3° ^ 
IH + 



a 



6« 



X 



a'- 



for X > 1 



for X < 1 



(4.3) 
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1 + ^ + O (^) for X = 1. To leading order in l/a, the solution of Eq. (| 



with ^ 

for X > 1 as a function of the conformal time has the following form: 

2 

1 I V \ 

a 



1 + Ztr 



JL 



(4.4) 



where r]tr is the conformal time corresponding to the transition redshift ztr- For x < 1, let 
r/o be the present epoch conformal time and normalise the scale factor as a(r]Q) = 1; the 
corresponding solution of Eq. ( |4.3| ) is: 

1 



(4.5) 



Joining solutions ( [1.4| ) and ( ^.5[ ) in a{ritj-), we can link the transition and the present epoch 
conformal time to the transition redshift as follows: 

Vo - f?tr = Ztr ■ (4.6) 

Note that the relevant contribution to the ISW effect comes only from solution (|4.5| ). In 
fact: i) the background solution (4.4) corresponds to a CDM dominated Universe and ii) 
from (^^) for x > 1 the speed of sound is exponentially vanishing for a — t- oo so that we 
can reasonably assume that Cg ~ 0. Therefore, if we substitute Eq. (4.4) and Cg = into 
Eq. ( |4.1| ) we obtain the same evolution equation for the gravitational potential as the one 
in a CDM dominated Universe, see [^]. In this instance, neglecting the decaying mode, 
<!>' = and no ISW effect is produced. 

Write then Eq. (4.1) combined with Eq. ( |4.5| ) and, from the leading order in x^" in 
(O) for X < 1, Cg = q: 



^„ 3 (1 + q) ^, 

Defining y = 1 + r]o — r], we recast Eq. ([ 

$ ^ 4- 

y 



3(1 + a) ,2 ' 



(1 + r/o - r/)2 
|) in the following form: 
■3(1 + 0; 



$ = . 



(4.7) 



+ k'^a 



^> = 



where the dot denotes derivation with respect to y. Equation (|4 
in terms of Bessel functions: 



y 



(4.8) 

can be solved exactly 
(4.9) 



where Ci and C2 are arbitrary integration constants. For large values of the order, the 
Bessel functions can be asymptotically expanded as follows |37]: 



eky/ay\ 
3a + 2 J 



3a/2+l 



^/3^TCX+'2^T 



1 



+ 



a 



a" 



and 



-+i 



{ky/ay) 



ekrs/ay\ 
3a + 2 / 



-3q/2-1 



37ra + 2tt 



1 + - + 

a 



a" 



(4.10) 



(4.11) 
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To leading order in 1/a, we plug the above asymptotic expansions into Eq. ( [4 .91) and find: 



3q/2 yZa 



2a 



vra 



-3«/2 



V3 



7ra 



(4.12) 



We assume the following initial conditions on the potential <I> in r/ = r/tr : ^ (vtr, k) = ^trik) 
and <!>' {rftr, k) = 0. The reason for this choice is that up to r/ = r/tr the potential behaves 
like in a CDM dominated Universe, i.e. it is constant. Solution ( [1.12 ) then becomes: 



1 



^tv{k) 1 - 3a \l + zt 



3a 



3a 



y 



1 — 3a 1 + ztr 



(4.13) 



For the calculation of the integrals ( |2.1| ) and ( p. 2D consider only the contribution propor- 
tional to y, i.e. 

^ M y (4 U) 

^tr{k) 1 + Ztr ' ^ ' ' 

since it is the dominant one for a — )• oo. Moreover, assume that the primordial power spec- 
trum Apj^ is the Harrison- Zel'dovich scale invariant one and that it propagates invariated 
up to 7?tr- For convenience, define 



D 



27r2 



25 



R > 



(4.15) 



combining the integrals and ( |2.2| ) with the approximated solution ( [4.14 ) and changing 
the integration variable from the conformal time to the redshift we write: 



1(1 + l)Cpw ^ 8;(; ^ ;l) 



27rL» 

Taking into account that ii{kz) 
lil + l)C7pw ^ 4^/(/ + 1) 



°° dfc 
T 



(4.16) 



211 J«+i/2(^-2)) we write Eq. ( 4.16 ) as follows: 



°° dA; 



J«+1/2(MJ«+1/2(M 



27rZ? (1 + z,,y 

Consider now the following case of the Weber-Schaflieitlin type integrals [37| 

r(/) 



(4.17) 



t2 



1 3 62 

4a'-V2r(« + 3/2)r(3/2) V 2' ^2'a2 



(4.18) 



where F is the Gauss hypergeometric function. Note that formula ( [4.18] ) holds true only 
if 5 < o. We perform the k integration in Eq. ( [4.17| ) according to Eq. ( [4.18[ ) and find the 
following expression: 



l{l + l)C7/sw _ + 1) 



27rL» 



(1 + ^tr) 



du , , , 

— Y / ClV V 



+ , (4.19) 

r(/-F-3/2) V 2' 2' u2 ' ' ^ ' 



where we have modified the integration range of v in order to satisfy the condition v < u and 
thus be allowed to apply Eq. ( 4.1^ ). Being the integrand function of Eq. (4.17) symmetric 
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with respect to the hne v = u, in Eq. ( [4. 191) we have recovered the correct value of the 
integral by multiplying by a factor 2. 

Expanding F in a hypergeometric series and performing the u and v integrations, we 
find the following series expansion for the ISW contribution: 

/(/ + l)Cpw _ ^ 1)^2 ^ r(/ + n)r(n - 1/2) 

2TrD (l + ztr)2 ^G + 2n + l)r(n + l)r(/ + n + 3/2) ' ^' ' 

which can be recast in the following more compact form: 

/(/ + i)cpw 20F4 r(/ + i) _1 [+1 ,^3 [+3 \ 

27rD ~ (l + ztrfni + m \ 2' 2 '' + 2' 2 'V' ^ 

where 3F2 is a generalised hypergeometric function |35]. 

The asymptotic behaviour of Eq. (|4.2l| ) for large values of I has the following form: 

+ ~ , (4.22) 

which is computed in Appendix It can be also directly obtained from the calculations 
of Kofman and Starobinsky (see Eq. (12) of [p6t ). 

In Fig. ^ we plot Eq. ( f4.21| ) and the relative asymptotic expansion Eq. ( 4.22| ) as 



functions of / for Ztr = 0.22 and ztr = 0.79. The former value is the best fit from the SNIa 
data analysis performed in the next section. 

The black line in the lower panel of Fig. |^ could be associated to those drawn in 
Fig. 5 of [25 1, where the ISW effect contribution is computed for the gCg up to a = 10. 



Note that the a — )• 00 contribution is smaller than the a = one, which corresponds to 
the ACDM case and was computed for the first time by Kofman and Starobinsky p6|| . 
We expected this result for the following reason. When a — )• 00 the Jeans wavenumber 
diverges. Therefore, according to ||2^, the contribution to Eq. (f4.16|) proportional to the 
fourth power of the speed of sound does not exist and the behaviour of l{l + l)C/^^/27r 
depends principally on the background evolution. The latter is pretty much similar to the 
ACDM one. The relevant difference is that for the gCg a — )• 00 model the transition from 
the CDM-like phase to the dS one is very sharp whereas for the ACDM is much smoother. 
Since in the former case the CDM-like phase lasts longer, the intensity of the ISW effect is 
lesser. 

Finally, note also how in Fig. 5 of |^] the trend of a decreasing ISW effect contribution 
can be already distinguished for a = 10 and large Vs. 

5. Background evolution for large a and SNIa data analysis 

In this section we address more quantitatively the behaviour of the background evolution 
for large values of a. Let us consider Friedmann equation ([4.3| ) to leading order in a: 

3 



o ~ < 



. (5.1) 

1 for Z < Ztr 
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Figure 2: Plot of the ISW contribution to the angular power spectrum l{l + 1)Ci/2-k normalised 
to D (solid lines) and of its asymptotic form for large Vs (dashed lines). Upper panel: ztr = 0.22. 
Lower panel: ztr ~ 0.79. 



As we pointed out in the previous section, Eq. ( p.l[ ) mimics the expansion of a pure CDM 
Universe for z > Ztr and the one of a dS Universe for z ^ Ztr- Note that, within this scenario, 
the present equation of state parameter is wq = —1, in contrast with the wq ~ —0.7 of the 
ACDM model. 

We now analyse the background evolution given by Eq. ( ^.1| ) on the basis of the 157 
nearby SNIa of the Constitution set |3S]. The supernovae data consist in an array of 
distance moduli defined as: 

^ = m-M = 51og(^^^ +25, (5.2) 

where m and M are, respectively, the apparent and the absolute magnitudes and Di^ is the 
luminosity distance expressed in Mpc units: 

, dz' c{l + z) r dz' , , 

DUz) = cil + z) =±-^ , 5.3 

Jo H{z') Ho Jo E{z') 
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Figure 3: Plot of the PDF vs the parameter ztr after margmalization over h. 



where E{z) is the Hubble parameter normahsed to Hq. 

The integral in Eq. ( |5.3D can be exactly solved for the Hubble parameter given in 



Eq. (|5l|): 



dz' 
E{^) 



' (3Ztr + 2) - 2 (1 + ^1 ^ ^^^^3/2 f^j, ^ ^ 



(5.4) 



for Z < Ztr 



Following Q, we assume flat priors for ztr and h and assume that the distance moduli are 
normally distributed. The probability density function (PDF) of the parameters has then 



the following form |40|: 



where [lo is the set of the observed distance moduli, 



(/i,^;tr) 



E 

i=l 



//o,i-51og (i?L/Mpc)-25 



and the normalisation constant C has the following form: 

1 



Azi 



■tr 



d/i e 



-x2(/l,2tr)/2 



(5.5) 



(5.6) 



(5.7) 



where the integration ranges over the parameters are, in principle, h € (— oo, cx)) and ztr £ 
(— l,oo). However, we choose the more reasonable ranges Ztr £ (0,2) and h G (0.5,0.9). 

In Eq. (5^) cr^^ , are the estimated errors in the individual distance moduli, including 
uncertainties in galaxy redshifts and also taking into account the dispersion in supernova 
redshifts due to peculiar velocities. 

After marginalization over /i, i.e. integrating Eq. (|5.5| ) in G (0.5,0.9), we show in 
Fig. ^ the PDF for the parameter ztr. The most probable value is ztr = 0.222. At the 68% 
confidence level ztr G (0.198,0.246), at the 95% confidence level ztr G (0.174,0.270) and at 
the 99% confidence level ztr G (0.154,0.290). 
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6. Summary and Conclusions 



In the present paper we have investigated the production of ISW effect within the gener- 
ahsed Chaplygin gas cosmological model. Thanks to an argument based on the Meszaros 
effect it is possible to find the new constraint a > 350. For this range of values, the Jeans 
wavenumber is sufficiently large so that the resulting ISW effect is not strong. Indeed, 
through a direct calculation, we have found a confirmation of the above qualitative con- 
straint because in the limit a — >• oo the ISW effect contribution to the CMB angular power 
spectrum is very similar to the one computed for a = 0, i.e. for the ACDM model. 

We have then addressed the background evolution of the Universe for a — t- oo and we 
have found that the model behaves like CDM at early times and then abruptly passes to 
a dS phase. Taking advantage of the SNIa Constitution set analysis, we have placed the 
transition at a redshift Ztt = 0.22. 

In conclusion, it seems that the gCg model has some chances of being viable not only 
for very small values of a but also for very large ones (we have here limited our discussion 
to the ISW effect only). However, it must be pointed out that in both cases a degeneracy 
problem appears: i) for a — )• 0, it is well-known that the gCg model degenerates into 
the ACDM; ii) for a — )• oo the degeneration takes place into a "step-transition" CDM- 
dS model. Note that in the second case, the degeneration is not complete. In fact the 
"original" CDM-dS model has a vanishing speed of sound. In the corresponding limit of 
the gCg, instead, the speed of sound diverges, so the scenario is completely different. A 
more complete analysis of such picture would therefore be interesting and could perhaps 
be performed in the framework of the Cuscuton model introduced in |41, 42]. 
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A. Calculation of the asymptotic behaviour of the ISW effect contribution 
for large /'s 



Consider Eq. ( f4.21| ), recast in the following way: 

^G + i)gp^ ^ 2V?4 rq + i) (l + l 1 3 / + 3 

27rD (l + z,^fTil + 3/2)''\ 2 2'^ 2' 2 ' 

Adopt the following transformation for the 3F2 function [^]: 



(A.l) 



3F2 (a,6,c;e,/;l) 



r(a)r(s + 6)r(s + c) 



3F2 {e - a, f - a, s; s + b, s + c;l) , (A.2) 



where s = e + f — a — b — c and {a, c, e, /} are arbitrary parameters. Employ Eq. ( [A.2p 
for the generalised hypergeometric function in Eq. (A.l) and find: 



i{i + i)cr _ 2^zl r(^)r(/ + i)r(3) / U,..^ 

2nD (1 + r (^) r(5/2)r(/ + 3) V ' 2 + 2'' + ' 



(A.3) 



Making explicit the series expansion of the generalised hypergeometric function we obtain 

^(/ + i)cpw_ 2V^4 r (^)r(/ + i) - r(i + i + n)r(n + 3) 



27rL> 



(1 + ztr)' r (^) r(i + 1) r(n + 5/2)r(/ + 3 + 



We approximate the Gamma functions for large Vs by means of Stirling formula [^] and 
find 



/(/ + l)Cpw 1 



27rL> 



V- 



r(n + 3) 



(1 + ztrf ^ ^o2"r(n + 5/2) 



(A.5) 



The series in Eq. ([A.5|) can be summed: 



t r(» + 3) 



n=0 



2"r(n + 5/2) 



and from Eq. (|A.5D we obtain the result of Eq. ( |4.22 ). 



(A.6) 
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